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S  U  M  M  A  R  Y 


A  steady-state  mathematical  model,  which  computes  the  spatial  configura¬ 
tion  and  tension  along  an  orbiting  towline,  has  ueen  developed  and  is 
described  and  presented  herein.  Numerical  results  are  in  agreement  with 
existing  flight  test  data.  Parametric  studies  of  the  operational  TACAMO  II 
system  indicate  that  very  high  verticaiity,  of  the  order  of  75  percent  or 
more,  can  always  be  achieved.  Numerical  results  presented  herein  Indicate 
that  significant  improvement  in  verticaiity,  time  on  station,  fuel  economy, 
and  flight  crew  comfort  can  be  achieved  in  the  TACAMO  systems  by  an  optimiza¬ 
tion  based  upon  the  mathematical  model. 

The  model  revoals  the  existence  of  multi-valued  .olutions  (alternate 
spatial  configurations)  under  certain  operating  conditions.  Sudden  tran¬ 
sition  from  one  steady-state  configuration  to  another  causea  high  transient 
dynamic  loads  .  Such  sudden  jumps  were  experienced-  experimentally  in  the 
early  flight  tests  of  TACAMO  II.  Operations  to  circumvent  these  potentially 
troublesome  regions  have  been  developed,  based  upon  numerical  results  of 
the  mathematical  model.  The  potential  for  achieving  increased  verticaiity 
by  usage  of  special  towlines  with  lower  dynamic  drag  coefficient  are  discussed. 
A  full  scale  lnotrumented,  flight  test  program  to  verify  the  steady-state 
model  and  the  development  of  a  transient  model  in  accordance  with  criteria 
presented  herein  are  recommended. 
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I.  INTRODUCTION 

Renewed  intereft  ha*  developed  in  the  dynamics  and  the  g-petlal  configura¬ 
tion  of  the  long  orbiting  tovllne.  Thin  problem  can  be  atudled  by  construct¬ 
ing  a  mathematical  model  which  practically  simulates  the  physical  situation. 
Studies  have  previously  bean  made  by  various  Investigators  (reference!  (a), 
(b),  (c),  and  (d)).  Each  study  yielded  Its  own  mathematical  modal.  Soma  ware 
overly  simplified  while  the  others  ware  overly  complicated.  Moreover, 
these  models  are  found  to  be  In  poor  agreement  with  experiment  or  not  yt 
verified  at  alii  a  satisfactory  mathematical  model  has  yet  to  be  developed. 

This  report  encompasses  a  critical  review  of  some  of  the  existing 
mathematical  models,  recomaended  criteria  for  mathematical  models  for  the 
long  orbiting  towline,  and  a  presentation  of  a  steady-state  mathematical 
model  developed  at  NAVAIRDEVCEN.  (Naval  Air  Development  Center)  under 
reference  (e) . 

II.  REVIEW  OF  PREVIOUS  MATHEMATICAL  MODELS 

A  literature  review  has  revealed  the  results  of  several  previous 
Investigations  concerned  with  the  mathematical  analysis  of  orbiting  tow- 
lines  (references  (a)  (b)  (c)  and  (d)).  Each  previous  Investigator  developed 
a  model.  Generally,  these  models  can  ba  classified  In  two  categories,  the 
steady-state  and  the  transient,  k  steady-state  model  simulates  a  towHne 
under  steady-state  motion.  Under  this  condition,  the  dynamic  and  geometric 
parameters  of  the  towline  are  Independent  of  time  except  for  either  a 
constant  speed  translation  at  a  constant  altitude  or  a  rigid  body  rotation 
about  an  axis  perpendicular  to  the  earth's  surface.  Mathematical  nude  Is 
to  simulate  situations  other  than  steady-state  motion  will  be  referred  as 
transient.  These  models  have  been  reviewed  and  are  discussed  Individually 
below. 


A.  Douglas  Aircraft  Company  published  a  report  In  February  1957 
titled,  "Some  Calculations  on  the  Long  Line  Technique  for  Lowering  Sonar 
from  Flying  Aircraft"  (reference  (a)).  It  is  a  steady-state  orbiting  towline 
model.  Even  though  the  purpose  of  that  work  was  different  from  the  current 
objective,  the  basic  dynamics  and  the  Interest  In  equilibrium  configurations 
were  the  same.  The  aerodynamic  force  acting  on  the  line  Is  assumed  to 
follow  the  sine  square  law  for  an  Infinitely  long  cycllnder.  Skin  friction 
In  the  direction  of  the  relative  velocity  was  not  considered.  The  solution 
for  the  equilibrium  configuration  of  the  cable  was  obtained  through  numerical 
integration  of  the  force  equilibrium  equations.  Integration  was  commenced 
from  the  tow  object  end.  Since  the  mcdel  Is  a  steady  ntate  model,  It  requires 
these  assumptions: 

1.  No  wind 

2.  Towplane  flies  a  perfect  circular  pattern  at  a  constant 
attitude,  and 

3.  No  cable  reel-in  or  reel-out. 
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in  the  formulation  of  the  differential  equation  it  ia  assumed  that  the  cable 
la  inextenalble  but  perfectly  flexible  in  bending. 

B.  In  1962,  Grunman  Aircraft  Engineering  Corporation  (reference  (h)) 
developed  a  mathematical  model  to  determine  the  dynamic  configuration  of 
a  long  tovllne.  Again,  the  purpose  of  that  model  was  different  from  the 
current  interest.  It  was  developed  to  study  a  long  line  orbit  technique 
for  surface-to-air  retrieval.  This  model  can  be  used  to  acquire  transient 
towline  configurations  Including  tha  effects  of  wind.  The  following  assump¬ 
tions  were  made  in  the  mathematical  formulation: 

1.  Aerodynamic  forces  acting  on  the  cable  follow  the  sine  square 
law  with  no  skin  friction  present. 

2.  The  towline  deforms  elastically  under  tension  but  it  cannot 
be  subjected  to  compression.  Also,  the  towline  is  perfectly  flexible  In 
bending. 


3.  At  the  upper  end  of  the  towline,  motion  is  specified  and 
at  the  lower  end  a  point  mass  Is  assumed. 

Grumman' s  modal  takes  account  of  the  longitudinal  elasticity  of  the  towline. 
Because  the  longitudinal  elastic  wave  speed  is  high  (compared  with  the 
lateral  wave  propagating  along  the  cable)  very  small  time  Increments 
(0.001  sec.)  are  required  between  successive  calculations.  This  makes 
the  computation  time  required  to  complete  a  realistic  initial  value 
problem  prohibitively  long,  even  with  the  aid  of  the  electronic  computer. 
Therefore,  primary  results  presented  by  Grumman  were  obtained  through  a 
simplified  version  of  the  modal  assuming  that  tha  towline  does  not  stretch 
with  tension;  It  Is  longitudinally  rigid.  The  inclusion  of  longitudinal 
flexibility  Is  unwarranted  unless  the  cable  is  subjected  to  rapid  varia¬ 
tion  of  longitudinal  loading  such  as  Impact  loads.  For  usual  applications 
the  complexity  is  not  justified  by  added  accuracy.  Tha  mathematical  tech¬ 
nique  used  In  constructing  the  physical  problem  is  the  lumped  mass  method. 

In  this  analysis,  the  cable  lc  divided  into  a  finite  number  of  segments. 

At  the  lower  end  of  each  segment,  a  concentrated  mass  and  a  projected  area 
equivalent  to  that  of  the  segment  is  attached.  These  point  masses  are 
assumed  to  be  connected  with  massless  straight  lines.  The  number  of  the 
straight  line  segments,  N,  which  can  be  used  to  approximate  a  continuous 
towline  Is  severely  limited  since  the  number  of  simultaneous  equations  to 
be  solved  Is  proportioned  to  N.  The  maximum  number  of  segments  used  in 
Grumman's  calculations  was  fourteen.  This  model  is  considered  to  be  in¬ 
adequate  for  the  purpose  of  the  current  study  where  accurate  towllne  config¬ 
uration  is  required. 

C.  Another  very  general,  and  complicated,  mathematical  model  of  the 
towllne  was  reported  by  the  Hayes  International  Corporation  in  1967 
(reference  (c)).  Hayes  studied  the  complete  maneuvering  tow  system,  which 
Includes  the  tow  plane  and  a  target  as  well  as  the  towllne.  The  mathematical 
technique  used  is  the  method  of  characteristics.  In  general,  this  mathe¬ 
matical  technique  yields  very  accurate  numerical  results  for  problems  of 
the  current  type  which  Involve  hyperbolic  differential  equations.  The 
accuracy  of  the  numerical  solution  could  have  been  further  enhanced  if 
center  difference  equations  were  used  Instead  of  the  forward  difference 
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equations  used  In  this  model.  As  In  most  of  the  oth#1*  models,  the  towline 
Is  assumed  to  be  inextenslble,  but  perfectly  flexible  In  bending.  Solutions 
of  the  equations  of  motion  of  the  towplene  and  the  target  were  used  as 
boundary  conditions.  These  boundary  conditions  Involved  fourteen  variables 
at  the  target  end  and  eleven  variables  ut  the  towplene  end.  It  Is  believed 
that  such  an  approach,  although  theoretically  more  accurate,  is  unwarranted 
in  this  problem.  The  added  accuracy  does  not  justify  the  greatly  Increased 
complexity  of  numerical  computations.  Computation  would  be  very  much  easier 
if  simplified  boundary  conditions  were  specified.  Since  the  towline  exerts 
only  a  secondary  effect  on  the  towplene,  the  motion  of  the  aircraft  may  be 
specified  at  one  end,  and  an  unpowered  drogue  of  simple  configuration  can 
be  used  at  the  other  end.  Such  a  model  would  suffice  for  the  study  of  the 
steady-state  long  orbiting  towline.  TVo  conditions  which  are  considered  to 
be  essential  In  the  study  of  transient  or  unsteady  towline  dynamics  were  not 
Included  In  Hayes'  mathematical  model.  They  are: 

1.  Wind  and  wind  shear 

2.  Towline  reel  rate 

Apparently,  Hayes  encountered  numerous  difficulties  it,  making  their  computer 
program  operational.  Included  in  their  report  were  numerical  results  from 
a  debugging  run.  It  is  believed  that  the  model  Is  not  yet  operational.  It 
was  not  at  the  time  reference  (c)  was  prepared. 

D.  As  part  of  a  study  on  an  orbiting  antenna  for  airborne  VLF  communi¬ 
cations,  Cornell  Aeronautical  Laboratory  (reference  (d))  developed  a  model 
to  analyze  the  motion  of  an  orbiting  antenna  cable.  This  model  1.  also  a 
steady-state  model.  The  aerodynamic  force  Is  assumed  to  follow  the  cross 
flow  principle,  but  the  skin  friction  was  assumed  In  the  direction  of  the 
cable  rather  than  in  the  direction  of  the  relative  air  flow.  This  assump¬ 
tion  Is  not  In  accordance  with  the  cross  flow  principle  as  defined  by 
Hoerner  (reference  (f)).  It  was  also  assumed  that  the  cable  Is  inextenslble 
but  perfectly  flexible.  The  equilibrium  configuration  of  the  antenna  Is 
established  through  integrating  the  dynamic  equilibrium  equation  from  the 
lower  end  (drogue  end).  The  dynamic  equilibrium  equations  are  written  in 
cylindrical  coordinates.  In  their  numerical  computations,  Cornell  neglected 
skin  friction.  Their  results  are  substantially  different  from  results 
obtained  by  the  NAVAIRDEVCEN  mathematical  model,  which  will  be  described 
In  detail  later,  particularly  In  the  region  where  multi-valued  solutions 
might  exist. 

III.  CRITERIA  FOR  MATHEMATICAL  MODELS 


After  critically  reviewing  these  mathematical  models,  it  can  be 
concluded  that  none  of  the  existing  models  can  be  used  for  the  current 
study  of  the  long  orbiting  towllnes  because  they  are  either  too  simple 
or  too  complicated  In  some  respect,  which  is  of  no  consequence  to  the 
present  Interest.  Therefore,  the  development  of  a  new  mathematical  model, 
which  adequately  simulates  the  dynamics  and  configuration  of  the  orbiting 
towline  under  various  operating  conditions,  was  undertaken.  The  method  of 
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characteristics  is  prsfsrrsd  to  the  lumped  mass  ms t hod  baeauss  It  generally 
givas  nova  aeeurata  numerical  rasulta. 

Practical  mathematical  models  for  th«  study  of  tha  dynamics  and  con¬ 
figuration  of  tha  orbiting  towline  should  satisfy  tha  following  criteria: 

A.  Criteria  for  tha  Steady- State  Modal: 

1.  Xt  should  be  able  to  yield  solutions  giving  tension  and 
geometric  configuration  along  the  entire  tovllna  for  the  following  operating 
conditions: 


a.  The  towplsne  travels  a  straight  and  level  path 

b.  The  tow  aircraft  travels  In  a  perfectly  circular  orbit 
at  a  constant  altitude. 

2.  Acceptable  Assumptions: 

a.  Aerodynamic  forces  follow  cross  flow  principle  (skin 
friction  In  the  direction  of  the  relative  air  flow) 

b.  Tovline  Is  Inextenslble  but  perfectly  flexible 

c.  No  wind  and  no  wind  shear  present 

d.  No  towline  reel -in  or  reel-out  motion 

e.  Air  density  varies  with  altitude  as  In  a  standard  day. 

3.  Boundary  Condition: 

a.  Towplane  flies  a  straight  line  or  a  perfect  circle  at 
constant  altitude,  flth  no  towline  being  reeled  In  or  out 

b.  Tow  drogue  is  of  simple  nose  or  center  of  gravity  tow 
configuration;  and  It  Is  not  powered. 

4.  Formulation  Should: 

a.  Admit  numerical  techniques  which  are  of  high  degree  of 

accuracy. 


b.  Allow  a  computer  program  that  can  provide  solutions  for 
many  different  operating  conditions  as  well  as  for  many  different  towline 
systems  In  one  computer  run. 

B.  Criteria  for  the  Transient  Model: 

1.  It  should  be  able  to  yield  time  dependent  solutions  for 
tension  and  geometric  configuration  along  entire  towline  for  the  follow¬ 
ing  operating  conditions: 
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a.  Th«  towplane  antari  Into  a  circular  orbit  Crow  •  toady* 
stata  straight  and  level  flight  with  the  towline  taing  reeled  out. 

b.  The  motions  of  the  towplane  changes  fron  s::e«iv  circular 
orbit  to  straight  and  level  with  the  towline  being  reeled  in. 

c.  The  towplane  undergoes  a  coordinated  turn  with  »  ct  rjitanc 
bank  angle  In  an  environment  with  non-uniform  s'  -'  v 1  nd . 

d.  The  towplane  undergoes  a  coord,  s'  d  ton  with  i-.:  risdy, 
non-uniform  wind  to  Indicate  stability  ofwulti-vs  4  editions. 

2..  Acceptable  Assumptions: 

a.  Aerodynamic  forces  follow  cross  flow  ^rinciplw  (s'.:  it. 
friction  in  the  direction  of  the  relative  air  flow). 

b.  Towline  is  inextensible  but  perfectly  flexible 

3.  Boundary  Conditions: 

a.  Towplane  velocity  and  reel  speed  are  prescribe! 
functions  of  time  with  very  low  rate  of  change  of  alrcri if  accei  ii  ’lions 

b.  Tow  drogue  is  of  simple  nose  or  center  o'  *r»iv..r\  tov 
configuration,  and  it  is  not  powered. 

4.  Formulation  should: 

a.  Admit  numerical  techniques  which  are  innr  tnti  o<  r. urate, 

thus  permiting  the  initial  value  problem  to  be  carried  out  repe.-.:d  y  for 
many  time  Intervals;  the  real  time  involved  a  transient  .  trie::  -?  b< 

of  the  order  of  minutes. 

b.  Allow  solution  with  a  computer  program  that  can 
executed  in  a  medium  slse  computer  with  reasonable  expenditures  ' 
computer  time. 

5.  Computer  programming  should  have  the  following  f  »xlt 

a.  Variable  towline  physical  and  aerodynamic  pr^ert1* 

b.  Variable  air  density  as  a  function  of  density  altit’*-‘*v 

c.  Variable  towline  length  between  data  points  fo~  bet l  > 
choice  of  time  increment  between  successive  calculations  as  well  as 
improving  the  accuracy  of  the  solution. 

d.  Redistribution,  addition  or  elimination  of  data  points 
along  the  towiine 
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IV.  NAVAIRDEVCBN  STEADY- STATE  MATHEMATICAL  MODEL 

In  the  study  of  the  transient  notion  of  e  towline  the  initial  configura¬ 
tion  oust  be  given.  Usually  the  notion  of  the  towlins  system  commence t  with 
a  constant  speed  straight  and  level  translation  or  a  constant  angular  apeed 
rotation  about  an  Imaginary  axis  perpendicular  tc  the  surface  of  the  earth. 
Sines  the  constant  speed  straight  and  level  translation  is  equivalent  to 
a  rotation  with  conn  tan  t  tangential  velocity  about  an  axis  at  infinity, 
a  simple  teadystate  model  which  describes  the  steady-stcte  notion  of  the 
tovline  orbiting  around  a  vertical  axis  will  provide  the  initial  data  for 
any  subsequent  transient  solutions.  On  this  premise,  a  steady-state  mathe¬ 
matical  model  which  satisfied  the  current  Interest,  and  which  will  be  referred 
to  as  the  NAVAIROKVCEN  model,  has  been  developed  and  is  described  below. 

A.  Assumptions 

Assumptions  used  in  constructing  the  mathematical  model  Include: 

1.  Towline  is  lnextenslble  but  perfectly  flexible 

2.  Aerodynamic  forces  acting  on  the  towline  follow  the  cross -flow 
principle  as  described  in  reference  (f) 

3.  Towplane  travels  in  a  perfectly  circular  path  at  a  constant 
altitude  with  no  towline  pay-out  or  reel -in. 

4.  Mo  wind  and  no  wind  shear  present 

5.  Air  density  varies  with  altitude  as  in  a  standard  day. 

B.  Formulation  of  Beuatlona 

It  is  convenient  to  formulate  the  problem  in  cylindrical  coordinates. 
Figure  1  shows  the  towplane,  the  towline  and  the  drogue  in  connection  with 
the  associated  coordinate  system.  The  equation  which  governs  the  steady- 
state  conflguretlon  of  the  towline  can  be  derived  as  follows  by  summing  up 
the  forces,  including  the  Inertia  force,  acting  on  an  element  of  the  tcwllneAS 
and  equatlng_£hc  sum  to  sero: 

jf*J  V-  fc(MroJjr*nuc  &>rce)4S  +  6r*\ri fatten 

force)**  +  fj(  Inert?*  force)* aj  sso  ( / ) 


The  aerodynamic  force  is  expressed  as  follows; 


where : 


CD  ■  Drag  coefficient 
^  ”  Air  density 
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Relative  air  velocity  vector 


«  Component  of  the  relative  air  velocity  In  the  direction 
perpendicular  to  the  towline 

d  «  Diameter  of  the  tovllne 

The  first  term  In  equation  (2)  Is  the  dynamic  drag  acting  on  a  unit  length 
of  the  towline.  The  direction  of  the  force  la  perpendicular  to  the  towllne 
with  a  magnitude  of  /•  */» 


The  second  term 


represents  the  skin  friction  which  Is  oriented  In  the  direction  of  the  rela- 
tlve  air  velocity.  Since  it^is  assumed  that  there  is  no  wind 

r?  _  _  (s,  t) 

Vfe/  -  TT~~ 

Referring  to  Figure  l,  it  is  clear  that  the  following  relationships  exist: 

r<s,t)  m  fi(s)  eR(Sjt)  +  £ 


+cojei 

(s,r)  s  cos  6  *  -  sih  £  J, 

where:  -  Q(S,0)+ult  *  &&)  + 


(u)ls  the  constant  speed  rigid  body  rotation) 


therefore, 


\n~~ 


Mb*  . 

TF  "^TF —  =  ^z-HO)uie9 


The  magnitude  of  *Rz!  Is: 

/&>,/  =  -MW 


Let  eg  denote  a  unit  vector  pointing  In  the  positive  direction  o£  ds. 

Ihe"' 

=  «.  <  0  4 * 

=  +  HU)  8'%  +  i'd)  -I 

where  prime  Is  used  to  denote  the  first  derivative  of  the  corresponding 
dependent  variable  with  respect  to  s. 

With  this  relationship,  the  component  of  the  relative  air  velocity  perpendicu¬ 
lar  to  the  relative  wind  can  be  written  as: 

=  *«)»>*«)  &)[*'«)  %  ) 
=  ro)  ufao  z'cO  60)  'it+C  *<s?e  '-/)$  ] 
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Hence,  the  magnitude  of  ^ls: 

!?J  -  ta) 


Since 


iU)l=  Mtfe'a? 

ConeequentXy: 

l%l  ■  ru> 

Substituting  ths  expressions  for  I fyelj  /^/  and  ^ 

In  Equation  (2)  yields: 

•7?  ■ »  <i  J  f  ✓  T^kW* 

+K6'e'-£]-irf^-{R\Dx%  (s) 

The  gravitational  force  vector  la  simply: 

^  3  ~/*3  *  w 

where yU,"  mass  per  unit  length  of  the  towline 
g  -  gravitational  acceleration 


The  Inertia  force  vector  Is: 

-  yU/SuI1^  (SO 

Flnslly,  the  expression  esn  be  resolved  into  three  oomponente. 

#  “  -  £ [M&  +m**% + 171**1 

«■«  /?/-  r 

(rK')'^rg'o%  +(r*6'yZ+  rw/ffy  t(ri')7 
r*6 A]£  *fi >t*6 ‘)* rt'siJQ-tbi-yZ  ((,) 


Substituting  Equations  (3),  (4),  (5)  and  (6)  In  Equation  (l)  than 
dividing  by  AS.  results  In  the  following  thraa  scalar  aquations: 

(rtz*)'-  WK  cb  V  '-**0*  6  (7) 

(rue’)  *  TK'e-vCi (ej 

(Ti’)'+  £ Vi 6 'iV  '-*V*  -o  (?) 

These  ordinary  differential  equations  can  be  solved  by  the  finite 
difference  method,  and  the  complete  solution  of  the  towline  can  be 
obtained  step  by  step  starting  from  the  drogue  end. 

The  appropriate  difference  equations  written  for  two  points,  1  and 
2,  a  distance  of  AS  apart,  can  be  put  in  the_form: 

S'-  Ct>  ff* •('•**&*  1? COL)aJ  O  *) 

rj^A~rM=[~  ( iO 

rz  2/-  T,  B,  '=  (-  Cc  I  f eJ/tfvV  )*S  (/  2) 

Subscripts  1  and  2  are  used  to  Indicate  the  corresponding  values  at 
point  1  and  2,  respectively.  Bars  placed  above  the  variables  represent  . 
the  arithmetic  average  of  the  values  of  these  variables  at  points  1  and 
2.  To  be  more  specific,  T^  Is  the  tension  at  point  1  and 

?-^L 

Equations  (10),  (11),  and  (12)  are  good  approximations  for  the  original 
differential  equations  provided  AJ  is  small  and  the  second  and  higher 
order  derivatives  of  the  dependent  variables  are  small  .  Similar  approxi¬ 
mations  give: 

=  /et 

-  i,  +  i'*s 

and 

$2-  &,+  6  C/S} 

At  point  2,  one  can  write: 

h,l=  /  Ob) 

If  all  the  variables  at  point  1  are  known,  then  the  unknown  variables 
at  rolnt  2  can  be  determined  from  equations  (10)  through  (16).  Since  the 
syr tern  of  algebraic  equations  Is  highly  non-linear,  probably  the  easiest 
way  to  solve  them  Is  by  numerical  iterations.  Straight  forward  iteration 
schemes  employed  can  be  found  In  Appendix  I. 


(/3) 

(/+) 
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The  convergence  of  the  Iteration  scheme  le  rather  fast.  In  general. 
It  converges  in  two  to  three  iterations.  Boundary  conditions  at  the 
drogue  end  can  be  expressed  as  : 

(tR'\  =  -  ■W  ufitft)  f  -  ^ 

(/»£.  « 
where : 

Mass  of  the  drogue 

Fg  “  Aerodynamic  side  force  acting  on  the  drogue 

Fq  -  Aerodynamic  drag  force  acting  on  the  drogue 

13  Aerodynamic  lift  force  acting  on  the  drogue 

The  aerodynamic  forces,  Pg ,  7q  and  F^,  acting  on  the  drogue  are 
functions  of  ,  the  angle  of  attack  of  the  drogue,  and  ,  the  slide 

slip  angle.  These  angles,  d  and  ,  can  be  determined  from  moment 

equations  which  require  that  under  equilibrium  the  summation  of  moments 
acting  on  the  drogue  about  any  axis  be  sero.  Aerodynamic  characteristics 
for  the  conical  drogue  used  In  Appendix  I  computer  runs  are  derived  from 
data  presented  in  reference  (g),  with  the  drag  coefficient  taken  to  be 
0.6. 


C.  Description  of  Computer  Program 

Numerical  computations  were  performed  In  a  CDC3200  computer. 

A  listing  of  the  computer  program  is  given  In  Appendix  I.  Primary  data 
Inputs  required  are: 

1.  Towplane:  Altitude,  absolute  speed  and  turn  radius. 

2.  Towline:  Mass  per  unit  length,  diameter,  aerodynamic  drag 
coefficients  and  length. 

3.  Drogue:  Weight,  base  area  and  aerodynamic  characteristics 

To  start  the  calculations,  it  is  necessary  to  assume  a  vertical 
separation  between  the  drogue  and  the  towplane  and  a  turn  radius  of  the 
drogue.  With  the  angular  speed,  uu,  given  by  the  ratio  of  aircraft  speed 
divided  by  Its  turn  radius,  calculation  can  be  continued  from  the  drogue 
end,  step  by  step  up  to  the  aircraft  end.  Most  likely  either  the  calcu¬ 
lated  towplane  turn  radius  or  the  calculated  altitude,  or  both,  are 
different  from  those  values  given  in  the  input  data,  i.e.,  the  values 
of  the  given  problem  to  be  solved.  The  differences  between  the  calcula¬ 
ted  towplane  altitude  and  turn  radius  and  their  corresponding  values 
given  by  the  problems  are  used  to  choose  new  values  for  drogue  altitude 


10 


NADC-AM-6849 


and  turn  radius  to  be  used  In  the  next  round  of  computations.  Hits  process 
of  adjusting  the  drogue  altitude  and  turn  radius  continues  until  both  the 
calculated  altitude  and  the  radius  are  very  close  to  their  corresponding 
values  given  by  the  problem. 

Finally,  numerical  results,  R,  2,  6,  T,  R' ,  2'  and  R'  are  printed  at 
selected  points  along  the  towline.  A  typical  computer  output  is  exhibited 
In  Appendix  II. 

D.  Numerical  Results  for  TACAMO  System 

Among  the  numerical  results,  those  of  particular  Interest  for  the 
TACAMO  system  are  the  maximum  tension  at  the  tow  point  (towplane  end  of  the 
towed  line)  and  the  vertlcallty  of  the  line,  l.e.,  the  ratio  of  vertical 
separation  distance  between  the  towplane  and  the  drogue,  expressed  as  a 
percentage  to  the  total  line  length.  Vertical  separations  versus  towplane 
turn  radius  for  the  TACAMO  system  under  various  operating  configurations 
are  plotted  in  Figures  2,  3,  and  4.  Each  curve  represents  vertical  separa- 
tlons  for  various  aircraft  turn  radius  while  all  the  other  parameters  are 
kept  constant. 

Figure  2  reveals  that  the  vertical  separations  as  well  as  the 
towline  tensions  at  the  tow  point  are  single  valued  for  a  specific  operating 
condition.  The  towplane  orbiting  velocity  Is  relatively  low  (V  ■  150  knots). 
As  the  towplane  orbiting  velocity  Increases,  these  values,  vertical  separa¬ 
tion  and  tension,  tnay  become  multi-valued.  Indicating  that  It  is  possible 
for  the  towline  to  assume  different  steady-state  configurations  with  the 
towplane  flying  at  exactly  the  same  altitude,  speed  and  turn  radius.  Which 
steady-state  configuration  the  towline  will  assure  depends  on  the  stability 
of  the  various  configurations  and  the  transient  operating  conditions  of  the 
towllne  system  prior  to  reaching  the  steady-state  cjndltion.  Three  important 
questions  are  of  interest  in  this  respect: 

1.  Are  these  steady-state  configurations  stable? 

2.  Under  what  conditions  will  the  towline  move  from  one  steady- 
state  configuration  to  another? 

3.  Will  it  jump  from  one  extreme  configuration  to  another  without 
assuming  the  Intermediate  equilibrium  configuration? 

These  questions  cannot  be  answered  by  the  steady-state  mathematical 
model.  A  properly  designed  transient  model,  however,  may  provide  these 
answers . 


The  existence  of  multi-valued  configurations  was  apparently  first 
suggested  by  Cornell  in  reference  (d).  Since  Cornell's  model  did  not  Include 
the  effect  of  skin  friction,  their  results  should  correspond  to  those  shown 
by  dotted  curves  in  Figure  2.  It  is  clear  that  these  dotted  curves  are  quite 
different  from  the  corresponding  solid  curves,  which  account  for  skin 
frictions;  particularly  In  the  regions  where  the  multi-valued  solutions 
are  anticipated.  Ihus,  neglect  of  skin  friction  extends  the  region  of 
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multi-valued  configuration!  and  increaaaa  the  dlffarencaa  batwaan  the 
alternate  aolutlona.  Dlffarencaa  in  vertical  aeparation  and  tovllne  tcnalon 
ara  greater  without  akin  friction. 

Figure  3  showa  the  aarne  type  of  curves  for  three  different  towline 
lengths.  The  general  ahape  of  theae  curves  la  similar  for  the  sane  velocity, 
a  smaller  turn  radius  la  required  for  a  towline  of  shorter  length. 

Under  the  assumed  standard  atmospheric  conditions  air  density 
decreases  as  the  altitude  Increases.  Higher  vertical  separation  la  there¬ 
fore  attained  when  the  aircraft  is  flown  at  a  higher  altitude  provided  other 
parameters  are  kept  the  same.  This  effect  la  shown  In  Figure  4. 

Some  Interesting  characteristics  of  the  multi-valued  solutions 
were  observed  when  the  space  coordinates,  R,  Z,  and  8,  of  the  towline  were 
plotted  along  the  entire  towline  In  Figure  5.  It  la  Interesting  to  observe 
that  from  configuration  I,  the  lowest  vertlcality,  to  configuration  II, 
intermediate  vertlcality,  the  vertical  separation  and  the  angle  6  are  Increased, 
while  the  radius  13  reduced  for  every  point  along  the  entire  towline.  However, 
from  configuration  II  to  configuration  III,  the  highest  vertlcality,  the 
vertical  separation  Is  further  Increased  and  the  radius  Is  further  reduced, 
but  the  angle  remains  essentially  constant.  Determination  of  the  exact  signi¬ 
ficance  of  these  characteristics  with  respect  to  stability  will  require  the 
transient  model. 

K.  Comparison  with  Flight  Measurements 

A  few  experimental  data  points  acquired  from  flight  test  in  1963, 
reference  (g),  were  marked  in  Figures  2  and  3.  They  are  In  good  agreement 
with  the  numerical  result  obtained  from  the  NAVAIRDEVCEN  mathematical  model. 

Recant  flight  tests  conducted  by  the  NAVAIRDEVCEN  have  given  additional 
evidence  that  the  towline  tension  at  the  tow  point,  as  calculated  by  the 
mathematical  model.  Is  essentially  correct.  Results  of  those  tests  are 
compared  with  calculated  data  in  Table  I.  Since  the  objective  of  the  recent 
flight  tests  was  not  to  verify  the  mathematical  model,  tension  measurements 
were  made  under  conditions  not  precisely  corresponding  to  the  ideal  steady- 
state  conditions  assuoied  In  the  model.  In  a  few  Instances,  a  range  of  tensions 
was  recorded.  Therefore  the  comparison  presented  in  Table  I  ghoul-"  be  evaluated 
qualitatively  rather  than  quantitatively. 

Although  the  vertical  separation  as  well  as  the  tow  point  tei^ion 
are  essentially  verified,  no  experimental  determination  of  the  towline  con¬ 
figuration  between  the  aircraft  and  the  drogue  has  been  made.  Therefore, 
a  well  designed  instrumented  and  controlled  test  program  to  accurately  verify 
the  steady-state  mathematical  model  should  be  conducted.  The  experimental 
verification  program  should  answer  questions  such  as: 

1.  Is  the  assumed  aerodynamic  force  on  the  towline  adequate? 

2.  Do  m.<ltl-valued  solutions  exist?  What  are  the  stability 
characteristics  of  the  multi-valued  solution? 


12 


NADC-AM-6849 


F.  OthT  Numerical  Results 

Towline  eye  tern*  with  parameters  other  than  thoaa  of  the  TACAHO 
syatem  were  studied  for  the  purpose  of  obtaining  info  neat  ion  and  knowledge 
for  the  development  of  future  syetema.  Numerous  combinations  were  studied. 

Some  results  of  particular  interest  are  exhibited  graph lcelly  In  Figures  6 
through  10. 

Figure  6  shove  the  Influence  of  the  magnitude  of  the  coefficient 
of  dynamic  drag,  Cp,  ou  vertical  separation  and  towllne  tension.  It  is 
clear  that  the  same  vertical ity  can  be  obtained  with  a  largei  towplane  turn 
radius  (lower  bank  angle)  if  the  towllne  hes  a  smaller  dynamic  drag  coefficient. 
A  lower  drag  towllne  also  makes  transition  from  low  verticallty  to  high  ver¬ 
tical  ity  smoother  and  more  gradual  when  the  aircraft  enters  the  orbit  at  a 
reduced  turn  radius.  The  less  severe  transition  minimizes  the  possibility 
of  a  sudden  change  in  the  equilibrium  altitude  of  the  drogue  which  might 
develop  high  dynamic  load  and  cause  towllne  breakage.  One  way  to  reduce 
the  coefficient  of  dynamic  drag  la  by  weaving  hairs  into  the  towllne. 

According  to  the  General  Tire  Corp.,  reference  (h),  such  hairs  reduce 
cable  drag  to  one  half  of  that  without  the  hairs. 

If  the  towplane  performance  parameters  (airspeed,  altitude,  and 
turn  radius,  etc.)  for  a  particular  system  have  been  optimized  for  maximum 
towllne  verticallty,  additional  Increase  in  verticallty  may  require  reduc¬ 
tion  of  towline  coefficient  of  dynamic  drag.  Since  tovllne  drag  coefficient 
has  a  larger  effect  in  the  lower,  higher  curvature,  portion  of  the  towllne, 
it  may  be  that  low  drag  devices  will  be  applied  only  in  that  portion  of  the 
towline. 


The  effect  of  coefficient  of  skin  friction  drag  on  the  vertical 
separation,  as  a  function  of  towplane  turn  radius,  is  exhibited  in  Figure  7. 
This  effect  is  very  prominent  in  the  region  where  verticallty  changes 
abruptly  with  turn  radius.  A  higher  coefficient  of  skin  friction  results 
in  higher  towllne  tension,  but  makes  a  smoother  variation  of  verticallty 
with  turn  radius.  Incidentally,  the  haired  cable  will  provide  higher  skin 
friction  than  the  one  without  the  hair.  The  higher  skin  friction  also  helps 
in  the  event  of  a  sudden  change  in  vertical  separation  as  the  aircraft  enters 
its  final  orbit. 

The  easiest  way  of  eliminating  the  possibility  at  a  sudden  change 
of  vertical  separation  for  the  current  TACAMO  system  is  to  increase  the 
drogue  weight.  This  effect  is  shown  in  Figure  8.  Increasing  the  drogue 
weight  from  100  pounds  to  300  pounds  effectively  removes  the  multi-valued 
region  for  the  particular  case  studied.  It  is  believed  that,  due  to  the 
Increase  in  drogue  mass,  the  magnitude  of  variation  of  drogue  vertical 
position  under  practical  operating  conditions  may  also  be  reduced. 

Another  way  of  eliminating  the  possibility  of  a  sudden  change 
of  vertical  separation  for  the  current  TACAMO  system  is  the  application 
of  a  stepped  diameter  towllne,  such  as  that  shown  in  Figure  9.  The  upper 
third  of  this  towllne  (0.210"  diameter)  (connected  to  the  towplane)  and 
the  lower  third  (0.160"  diameter)  are  made  of  existing  cables  used  in 
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current  systems.  The  middle  third  ia  a  new  edible  of  0.184"  diameter.  Calculated 
vertical  aeparatlone  and  towline  tana Iona  for  the  etepped  towline  and  the 
correaponding  reaulta  for  the  existing  0.210"  towllne  ere  shown  In  Figure  10. 

The  stepped  diameter  towline  offers  the  following  potential  advantages 
over  the  existing  0.210"  constant  diameter  towllne: 

1.  No  multi-valued  solution  exists  In  the  operating  region 

2.  Higher  vertical lty  In  tha  low  vertical  ity  region  and  no 
reduction  In  vertlcallty  in  the  high  vertieallty  region, 

3.  Lower  towllne  tension  at  the  tow  point  (2180  pounds  Instead 
of  2840  pounds  at  75  percent  vertlcallty).  Because  of  the  low  towllne 
trail  angle  at  the  towplane,  this  Is  a  640  pound  reduction  In  tow  drag. 

Changing  the  base  area  of  the  drogue  has  little  effect  on  the 
vertical  separation,  as  shown  In  Figure  11,  where  the  base  area  of  the 
drogue  Is  varied  from  3.68  to  7.07  square  feet. 

Trends  due  to  variations  of  towllne  diameter  or  towllne  weight 
per  unit  length,  can  be  observed  in  Figure  12.  If,  as  the  towllne  diameter 
is  Increased,  the  mass  per  unit  length  is  also  Increased  In  direct  proportion 
to  the  square  of  the  diameter,  higher  towllne  vertlcallty  Is  obtained  and  the 
region  of  high  vertlcallty  variation  Is  shifted  to  the  right  (largest  turn 
radius).  However,  the  opposite  Is  true  lfr  as  the  diameter  of  the  towllne 
Is  Increased,  the  mass  per  unit  length  Is  kept  constant. 

Tension  at  the  tow  point  Increase*  as  the  mass  per  unit  length 
of  the  towllne  Increases.  Tension  remains  essentially  constant,  however,  if 
mass  per  unit  length  of  towllne  Is  kept  constant  as  the  diameter  of  the 
towllne  varies. 

G.  Discussion  of  Steady-State  Solutions 

1.  Avoidance  of  Region  of  Multi-Valued  Solution 

The  NAVAIRDEVCEN  steady-stats  mathematical  model  predicts  that 
there  are  operating  regions  where  multi-valued  solutions  exist,  but  the 
regions  are  not  as  wide  and  the  differences  In  vertical  separation  for  the 
extreme  solutions  are  not  as  great  as  those  predicted  by  the  Cornell  model. 

Until  more  Information  concerning  the  existence  and  the  stability  of  alter¬ 
nate  configurations  can  be  acquired,  It  Is  desirable  to  avoid  operating  a 
towllne  system  In  these  regions.  This  objective  can  be  accomplished  by  the 
following  procedures: 

a.  Increase  the  weight  of  the  drogue.  For  example,  the  TACAMO 
drogue  should  be  Increased  from  100  to  300  pounds. 

b.  Before  entering  the  orbit,  keep  the  towplane  speed  as  low 
as  possible,  using  flaps  if  available.  Gradually  Increase  bank  angle 
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and  Increase  towplane  velocity  slightly  si  higher  bank  angles  art  approached. 
Execute  the  transition  from  straight  and  laval  flight  to  final  orbital  flight 
slowly  in  no  less  than  three  complete  orbits. 

2.  Optimisation  of  Flight  Operations 

For  a  given  operating  system  with  fixed  towline  length,  there 
are  many  combinations  of  towplane  operating  parameters,  velocity,  turn 
radius  and  altitude,  which  will  result  in  an  orbiting  line  operation.  Ju¬ 
dicious  selection  of  the  towplane  parameters  can  optimise  the  operetlon  for 
a  number  of  different  purposes,  e.  g.,  greatest  towline  vertleallty,  minimum 
variation  in  vertleallty,  best  fuel  economy  for  the  towplane,  maximum  comfort 
for  the  towplane  crew,  etc.  These  optlmisstlons  can  be  easily  and  economically 
studied  with  the  aid  of  the  ateadv-state  mathematical  model.  Bold  steps  can 
be  taken  and  novel  ideas  can  be  explored  without  danger  to  the  towplane  crew 
or  equipment.  In  addition,  with  the  high  speed  electronic  computer,  numerous 
cases  can  be  calculated  quickly  add  Inexpensively. 

Review  of  scattered  field  operating  data  on  existing  systems 
indicates  that  significant  potential  Improvements  can  be  affected  by 
optimization  studies  on  current  operational  systems  as  well  as  on  future 
systems. 


3.  Optimization  of  Towlines 

Foregoing  portions  of  this  report  have  pointed  out  the  effects 
of  dynamic  drag,  diameter  and  weight  per  unit  length  on  the  behavior  of 
towlines  in  orbital  flight.  Low  drag,  large  diameter  and  heavy  weight  would 
appear  to  be  desirable  objectives  in  all  cases,  but  they  do  not  necessarily 
lead  to  the  optimum  towline  for  a  given  operation.  In  the  case  of  relatively 
short  towlines,  or  in  situations  in  which  the  towplane  speed  and/or  turn 
radius  cannot  be  sufficiently  reduced,  low  drag  and  high  unit  weight  are 
significant  in  achieving  vertleallty.  On  the  other  hand,  with  long  towlines, 
high  unit  weight  may  be  an  Intolerable  burden  on  the  towplane.  There  are 
many  possible  variations.  The  special  towline  of  stepped  diameter  represents 
one  oi  the  variations. 

The  NAVAIRDETVCEN  steady-state  mathematical  model  provides  a 
ready  tool  for  investigating  the  significant  characteristics  of  towlines, 
optimizing  them  for  current  operations,  and  defining  profitable  avenues 
for  the  development  of  improved  towlines. 

V.  CONCLUSIONS 

A.  No  previous  mathematical  model  adequately  represents  and  predicts 
the  characteristics  required  to  design  and  optimize  orbiting  towline  systems. 

B.  TWo  models  are  actually  required:  a  steady  state  model,  and  a 
transient  model.  Criteria  for  both  are  presented  herein. 

C.  A  steady-state  model  has  been  developed  at  the  NAVAIRDEVCEN.  The 
analytical  formulation  and  a  computer  program  are  presented  and  discussod 
herein. 
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D.  Numerical  results  for  a  number  of  cases  corresponding  to  operational 
orbiting  towline  systems  are  presented  and  discussed.  Theao  results  are  In 
close  agreement  with  existing  flight  test  measurements. 

E.  The  steady-ataee  model  predicts  multi-valued  solutions  under  certain 
operating  conditions.  These  regions  of  operation  should  be  avoided  until  the 
existence  and  relative  stability  of  the  alternate  solutions  are  determined. 
These  regions  can  usually  be  avoided  in  actual  operation, 

V.  Stepped  diameter  towlines  can  eliminate  the  multi-valued  region 
and  reduce  the  tow  load.  This  lessens  the  possibility  of  towing  breakage 
and  Improves  towplane  endurance. 

G.  Vertlcallty  of  orbiting  towllnea  can  be  increased  by  special 
towlines  having  reduced  dynamic  drag  coefficients. 

VI.  RBCOmBNDATIONS 

A.  The  steady-state  model  should  be  utilised  to  optimise  the  operation 
of  current  orbiting  towline  systems  and  to  design  and  optimise  future 
systems.  Significant  Improvement  In  vertlcallty,  fuel  economy,  time  on 
station,  and  crew  comforts  could  be  made  by  an  optimisation  of  the  current 
TACAMO  system  using  the  mathematical  model. 

B.  A  full-scale  Instrumented  flight  test  program  should  be  designed 
and  conducted  to  experimentally  confirm  the  numerical  results,  Including 
the  results  from  the  stepped  towline,  predicted  by  the  steady-state  model. 

The  existence  and  stability  of  the  multi-valued  solution  should  also  be 
examined. 


C.  A  transient  mathematical  model  should  be  developed  in  accordance 
with  the  criteria  presented  herein. 
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